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Abstract 

We study global existence problems and asymptotic behavior of higher-dimensional 
inhomogeneous spacetimes with a compact Cauchy surface in the Einstein-Maxwell- 
dilaton (EMD) system. Spacelike r^~^-symmetry is assumed, where D > 4 is space- 
time dimension. The system of the evolution equations of the EMD equations in the 
areal time coordinate is reduced to a wave map system, and a global existence theorem 
for the system is shown. As a corollary of this theorem, a global existence theorem 
in the constant mean curvature time coordinate is obtained. Finally, for vacuum 
^ ' Einstein gravity in arbitrary dimension, we show existence theorems of asymptoti- 

' cally velocity-terms dominated singularities in the both cases which free functions are 

analytic and smooth. 

o 

Cn ■ PACS: 02.30. J,., 0A.20.Dw, 04.20. Ex, 98.80. 

O 

^ ■ 1 Introduction and summary 
o ■ 

<sJ '■ Global existence problems of fundamental field equations must be solved as the first step of the strong 

^ I cosmic censorship, which states that generic Cauchy data sets have maximal Cauchy developments which 

are locally inextendible as Lorentzian manifolds. It is important to consider the strong cosmic censor- 
ship in scope of unified theories such that superstring/M-theory. As is well known, such unified theories 
■ predict that the spacetime dimension exceeds four and it is expected that such extra dimensions must 

rS I not be different from ordinary dimensions and should be taken into account in the asymptotic regions 

' (e.g. near the singularities) of spacetimes. Although higher-dimensional model have a long history within 

superstring/M-theory, mathematically rigorous results for cosmological spacetimes are much less under- 
stood than for stationary ones (for examples, see [GSW]). In particular, there has been less exploration of 
global in time problems in higher-dimensional, cosmological spacetimes. Recently, spatially homogeneous 
cosmological models was analyzed [HS] . Then, we want to study an inhomogeneous cosmological case as 
the next step. 

In the cosmological context, it seems convenient to consider spacetimes which develop from smooth 
Cauchy data given on a compact, connected and orientable Cauchy hypersurface. In addition, we would 
like to investigate inhomogeneous cosmological spacetimes with dynamical degree of freedom of gravity. 
Then, we will consider globally hyperbolic spacetimes (A^c g), with A4d = -Md-i x R a smooth D- 
dimensional manifold (Z? > 4), g a Lorentzian metric, which develop from smooth Cauchy data invariant 
under an effective action of Gd-2 = U{1) x U{1) x ■ ■ ■ x U{1) = T^^^ on a compact {D — l)-dimcnsional 
spacelike manifold Md-i- The same symmetry is assumed for matter fields if they exist. The resulting 
system of field equations becomes one of 1-1-1 nonlinear partial differential equations (PDEs). From 
Mostert's theorem [MP] it is admitted that M.d-i/Gd-2 is a circle and Md-i is homeomorphic to 
gi X ^ tD-i^ Therefore, we wih suppose Md-i ~ T^-^. 

In four-dimensional vacuum spacetimes with the above setting, Moncrief and Isenberg have proved 
global existence theorems of the Einstein equations in areal and constant mean curvature time coordi- 
nates [MV, IM]. Recently, it has been generalized to non- vacuum case [AH, ARW, HO, NM02, NM03]. 
One purpose of the present paper is to extend the previous results to higher-dimensional spacetimes. 

Once it has been shown global existence theorems, we want to know asymptotic behavior of cosmo- 
logical spacetimes as the next step. That is, it should be analyzed nature of spacetime regions near 
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singularity (if incomplete) and near infinity (if complete). In this paper, we will focus to consider nature 
of singularities. It is conjectured that, near a cosmological singularity, a decoupling of spatial points 
occurs in the sense that evolution equations cease to be PDEs to simply become, at each spatial point 
asymptotically, ordinary differential equations with respect to time [BKL]. In other words, cosmological 
singularities are locally asymptotically velocity-terms dominated (AVTD ^) or mixmaster. Concerning 
this BKL conjecture, it has been shown that four-dimensional Gowdy symmetric spacetimes with or with- 
out stringy matter fields (with dilaton couplings) have AVTD singularities in general in the sense that the 
singular solutions depend on the maximal number of arbitrary functions [KR, NTM] . These results have 
been made no-symmetric and higher-dimensional generalization [DHRW]. That paper has established 
AVTD behavior for vacuum gravity in spacctimc dimension Z? > 10 and for the Einstcin-dilaton-matter 
system with dilaton couplings in spacetime dimension D > 2. Thus, the question whether the AVTD 
behavior exists or not still remains open for vacuum Einstein gravity in dimensions D G [5, 10] rigorously. 
Another purpose of the present paper is to show that there is AVTD behavior for vacuum Einstein gravity 
in arbitrary spacetime dimension. This result complements the BKL conjecture by combining together 
the previous works mentioned above. 

The above results concerning AVTD behavior arc of C"^ (analytic) category. Recently, it has been 
extended to C°° (smooth) category in the case of four-dimensional vacuum Gowdy spacetimes with any 
spatial topology [RA, SF] . We will generalize this to higher-dimensional spacetimes and tie the result on 
AVTD singularities together with the global existence results. 

Let us consider a truncated action of the bosonic supergravity theory (i.e. low energy effective 
superstring/M-theory), which contains only the gravitational (metric), the dilaton and p-form fields. It 
can be shown that this truncation is consistent, in the sense that the fields that are retained are not 
sources for the fields that are eliminated [SK]. In this paper, we will consider the case p = 1. Thus, the 
system becomes the Einstein-Maxwell-dilaton (EMD) system. The action for the D-dimensional EMD 
system is given by 



Sd = J d^x^/^d^ [-^R + 2{d(Pf + e-^'"^F^ 



(1) 



where ^Ris the Ricci scalar with respect to g, is the dilaton field, i^^^ = 9^^^ — 9^^^ is the Maxwell 
field strength and a is a coupling constant. Varying the action we have the following field equations: 



D 



R^^ = 2d^<j)d^(t) + e 



D-2 



(2) 



D 



(3) 

(4) 



where /i, A run from to _D — 1, ^i?,, 



d[t,F,x] = 0, (5) 
,y and arc the Ricci tensor and the d'Alcmbertian with 



respect to g^i^, respectively. Note that truncation of the Maxwell fields or of both the Maxwell and the 
dilaton fields is consistent, but trivializing only the dilaton field does not induce the Einstein-Maxwell 
(EM) system from the the EMD system. In this case, as we can see from the dilaton equation (3), the 
EM system with a constraint F"^ =Q would be obtained. 



2 G£)_3-invariant spacetimes: Reduction of a wave map system 
coupled to three-dimensional gravity 

We take Aio-i to be a principal fiber bundle with compact, two-dimensional base E, with fibers the 
orbits of Gd-3 = U{1) x ■ ■ ■ x U{1) = T^~^, and with a metric g on Md invariant under the right action 
of Gd-3- Therefore, A^d is a principal fiber bundle with base S x R and group Gd-3- 
^In this article, the word "AVTD behavior" is equivalent to "Kasner-like behavior". 
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Let ^/ = di (/ = 3, • • • , D — 1) be (D — 3) spacelike commuting Killing vector fields, such that 
£jj£f = and = 0- Note that the Killing vectors are tangent to Md-i, hence spacelike. Form the 

hypotheses, the metric g for the spacetime is [CBDWM] 

g = f-'^-f^^dx^'dx'' + fijQ'e-', — dx^ + w'„,dx"', (6) 

where f^^jmn is a metric on E x R, m, n, i,j = 0, 1, 2, / = det//j > 0. 7rnn, W^m and fij are depend 
only on the coordinates x™, that is a necessary and sufficient condition for the metric g to be invariant 
under the right action of Gd-3- 

We consider the Maxwell field A^. Assume that there is a gauge such that the U{1) gauge field 
respects the spacetime symmetry C^jA = 0. Then we have djA^ = 0. According to [IW] (see also [lU]), 
we can define the Maxwell field potentials and ^ from equations (5) and (4) as follows. 

F^i = d^^i, (7) 

^-2a4>pmn ^ V ."^"^/^ VI/ (8) 

V— det7 

where e^nM denotes the permutation symbol such that eoi2 = 1 and tmni = 0. The equation (7) satisfies 
a part of the Bianchi identity (5), while the equation (8) fulfills a part of the Maxwell equations (4). 
From (7) and (8) the remaining Maxwell equations and the Bianchi identity can be described through 
$7, 4", fij, and 7mn- Introduce torsion defined by 



= ///Jy^d^e^,A.7""7^"5™M/^„ (9) 

where e/^o...^D_i is the Levi-Civita tensor eo...£)_i = 1 and 7'"" is the inverse metric of Jmn- One can 
obtain the following evolution equations for the Maxwell field potentials and ^ from equations (4) 
and (5): 

= fV^fijVm^K - /-^V"*a;/„e2«^ + 2aV™.^V„$/, (10) 

and 

m = f-^W'^fWm^ + /"V^^/Wj^e-^"^ - 2aV™(^V„*, (11) 

where V and □ = V^Vm are the covariant differential operator and the d'Alembertian of "fmn, respec- 
tively. 

The equations for the dilaton field takes the form 

□ = -^(e-2«<?'/^Jv"$/V„$.7 - e2"*'/-iV'"1'V„«'). (12) 

Now, we will write down the Einstein equations. From the mixed (/, m)-components of the Einstein 
equations (2) we have 

d^coi, - d^uif, = 2(9,$79^* - a^$/5,*). (13) 

From this, we can define the twist potential Qi such that 

dQi + iJ/ =cc>/ + $/d*-*rf$/, (14) 

where Hi is a harmonic one- form for some given Riemannian metric on S which is compact [CJ, CBM96]. 
For simplicity, we will suppose 

Hi = 0. (15) 

Evolution equations for Qi can be found by taking a three-covariant divergence of equation (14) taking 
account definition of w/^ (9) 

□Q/ = /-'(V"/ + W™*)(V„^Q7 + W„$/-$jV„^f) 

+ ^{^"'fij + */V"*j)(V„Qk + ^V^^K - $kV„*) 

+ /-i$,V™/V„^M/-/''^M/V'"//jV„$^. (16) 
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Next, evolution equations for the scalar multiplet fu follow from (7, J)-components of the Einstein 
equations (2): 

- 2e-2«^V'"$7V„$j + -^^//7(e-2«^/^^V"$KV„$i - e^^^Z-^V^^V™*). (17) 

The eifective Einstein equations on the three-spacetime whose metric is 7 follow from (to, n) and 
(/, J)-components of the Einstein equations (2): 

''Rmn = lif-^VmfVnf + f"f''''ymflKynfjL) + e-^"^^/" V„$/ V„$ J + e^'^^/"! V^M/V^* 

+ ^/"'/"(V^Q/ + *V„$7 - $/V„*)(V„gj + *V„$J - $7V„*) + 2Wm<l>^n<l>, (18) 

where ''Rmn is the Ricci tensor for -fmn- 

Note that the Maxwell equations (10) and (11) can be rewritten by using Qj as follows: 

= /■'''V^'fijVm^K - /-'^"^(V^Q/ + *V„$7 - ^iVm^)e'""t' + 2aV"<^V„$7, (19) 

and 

□ * = /-lV"/V„^^' + /"V'"$/(V„^OJ + *V„$j - $jV„*)e-2'^'^ - 2aV'"<^V„*. (20) 
When fij, -ymn, Qi, ^i, ^ and (p are known on S x R, we can get two- forms on S x R, where 
W/,. = diW'j - djW'i = -/-V"V-det7eym7""(5nQj + * 9n*j - *j9„*). (21) 

We can deduce from one-forms on A^^i if and only if the following two conditions hold [CJ, 

CBM96]. One is a local condition which is the system of evolution equations. Another is a global 
condition. Since each U (l)-symmetric one-form is independent, in the case S compact, Gd-3 = T^~^, 
the global condition reads 

ni=^ I dQ' = ^ ( {m e Z). (22) 

In summary, the following action describing a wave map coupled to three-dimensional gravity is 
obtained as an effective action: 



= j d^x^ 



-det7 



(23) 



where ^i? is the Ricci scalar for ^mn- It is easy to show that the evolution equations, (12), (16), (17), 
(19), (20), and the three-dimensional Einstein equations (18) are obtain from this action (23). The only 
wave map consisting of scalar fields has dynamical degrees of freedom, but three-dimensional gravity does 
not have ones. 

3 Global existence theorem for T^~^-symmetric spacetimes in 
areal time coordinate 

We will consider Gd-2 = U{1) x ■ ■ ■ x U{1) = T^^^-invariant cosmological spacetimes. That is, we 
assume the existence of another spacelike Killing vector field ^2 = ^2 which commutes with the other 
Killing vectors [^2, = 0. The Maxwell and the dilaton fields are also assumed to be independent of the 
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coordinate x^. Thus, M d-i can be parametrized by and — 6 £ M d-i/Gd-2 ~ [0, 27r]|mod27r ~ S^i 
where I = 2, •••,£)— 1. The Frobenius integrabihty condition is [DFC] 

=0. (24) 

Spacetimes satisfying the above condition admit a foliation by two-dimensional integrable surfaces or- 
thogonal to the Killing fields ^i. From equation (9) one can obtain that the condition (24) for I = / 
is 

w/2 = 0. (25) 

The equations (25) arc satisfied if they hold somewhere at one point in the spacetimes. Indeed, we have 
the following equation by (13): 

d^uji2 = 2{d2^id^^ - di^^A^) = 0, (26) 

because of 82^1 = 82^ = 0. Since it need not to distinguish ^2 from the Frobenius integrabihty 
condition (24) is satisfied for I if equations (25) hold at least one point ^. Hereafter, we assume this 
hypersurface-orthogonality condition. 

Under the above assumption, the line element of the spacetimes takes the 2 x 2-{D — 2) x (Z) — 2) 
block-diagonal form. Thus, the three-dimensional Lorentzian metric can be written in the form 

lij = e^^i-dt'' + de^) + p'^di;^, = Wl = 0, (27) 

where t = x^, ip = x"^, A = \{t,0) and p = p{t,0). Here, we have used the fact that two-dimensional 
spacetime (which is spanned by t and (^-coordinates) can be conformal flat. In this coordinate, we have 
a linear wave equation for p from the Einstein equations (18) for three-spacetime, 

{d^ - dl)p = 0. (28) 

By using equation (28) and the same argument with Gowdy [GR], we can take areal time coordinate 
p = t. Then, the spacetime metric is 

g = f-\e^^{-dt'^ + de^) + + fij{dx^ + d'ilj){dx'^ + W'^dip), (29) 

where the metric functions depend only on t and 9. The same dependence is assumed for functions of 
matter fields. Note that a metric of T^-Gowdy symmetric spacetimes would be induced from the metric 
(29) if £» = 4 [GR]. 

In the areal time coordinate (29), we have constraint equations from (18) as follows: 

2t-'dtX = Ir'mff + {dom + lf"f""'idc.fiMdtfjN + defiMdefjN) 



e 



-2a4> fIJ 



+ \f~'f" li^tQT + - ^idt^) [dtQj + - ^jdt^) 

+ [deQl + -^dg^i - ^idg^) {dgQj + <idg^j - ^jdg^)] 

+ 2[{dtct,f + {dgct>fl (30) 
+ e-^'"^f"dt^idg^j + e'^'"l'f-^dt^dg-^ 

+ \f~^f" i^tQl + ^dt^i - ^idt^) {dgQj + ^dg^j - ^jdg^) 

+ 2dt<f>dg(P. (31) 



^In the case of four-dimensional U{1) x {7(l)-symmetric spacetimes, it is known that the above assumption corresponds 
to vanishing twist constants in the case of T^-spatial topology. In addition, if there is a symmetry axis in the spatial section, 
lilfe 5^ X 5^, the condition (25) holds at the axis, then the hypersurface-orthogonality is automatically satisfied [CP]. 
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The integrability condition dtdgX = dgdtX of equations (30) and (31) is assured whenever the other 
evolution equations are satisfied. Note that we can obtain a evolution equation for A from equations for 
others. 

d^X - dlX = i/-2[_(a,/)2 + (a,/)2] + ^/"f'^'i-dtfiMdtfjN + defiMdefjN) 
+ {deQi + <^de^i - -^ide^) {dgQj + ^dg^j - ^jdg^)] 

+ 2[-{dt^f + {dgct>n (32) 

Thanks to the arcal time coordinate, the metric function A is decoupled with other fields. Therefore, it is 
enough to solve the evolution equations at first. After that, we must demand compatibility conditions for 
A. That is, equations (30) and (31) are ones to determine the metric function A. Under the coordinate 
(29), we will call a system of equations (12), (16), (17), (19), (20), (30), (31) T^-'^ -symmetric EMD 
system. 

Now, we have the following conclusion by direct calculation: 

Lemma 1 Let M = R x and N = Ri(''^-2)(D+i) manifolds with Lorentzian metric r] and Rieman- 

nian metric h, respectively. Then, a wave map U : M ^ N is equivalent to the evolution equations (12), 
(16), (17), (19), (20), of the T^~^ -symmetric EMD system. Here, the action for the wave map is 

SwM = j dtdOdtpC^M = j dtd9dtp^/-detr]r]°'^hABdaU'^d0U'^ , (33) 



and the metrics are 



and 



-dt^ + de"^ + t^dTp"^, < 6», V < 27r, (34) 



h := -^df^ + ^f^f'^'^dfiKdfjL + e-^'''^f"d^id<i>j + ^d^^ 

+ hf"{dQi + - ^id^){dQj + - $jd*) + 2#2. (35) 

2/ 

Note that U is independent ofip. □ 

The integrals (22) are independent of time t. This fact follows from Noether's theorem [BCM, SS]. 
Let Z he a Killing vector for the metric h given by (35). It is well known that the following quantity is 
independent of the choice of compact Cauchy hypersurfaces, 

E{Z,S)= j^rj^^hAB-^Z^dSp, (36) 

where q:, 3 = t, 9, ijj, = (P/, Q/, ^f, 0)^ and dSa = da\/ {dtAdO Ad4>), and V denotes contraction. It 
is easy to sec that arc Killing vectors for the metric h. Taking Z = and S ~ {t — constant} = S 
one obtains the integrals (22). Thus, they are conservation quantities when the wave map system holds. 

Although Lemma 1 is quite general and useful, geometry of the target space {N, h) is very complicated 
since scalar multiplet fu does not be fixed. Therefore, to analyze the system furthermore, we assume 
the following, 

fu = e^'^'Sij, (37) 

and P = Yli^i- This assumption docs not restrict to diagonal (i.e. polarized) spacetimcs. Indeed, our 
T-^~^-symmetric spacetimes with the assumption include four-dimensional unpolarized Gowdy symmetric 
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ones as a special case. Under the assumption (37), the metric of the target space can be written by 

7=3 7=3 

D-1 

+ -e-^-f" ^ e-"^^' {dQi + *di>/ - ^id^f + (38) 

7=3 

For this target space, we can show the following global existence theorem: 

Theorem 1 Let {AiD, g^^iy, A^^, (p) be the maximal Cauchy development of smooth T^^^ -symmetric Cauchy 
data on Md-i ~ T^~^ for the T^~'^ -symmetric EMD system. Then, under the assumption (37), A4d 
can be covered by the areal time coordinate with t £ (0, oo). 

It is known the local existence result for the wave map system (e.g. Theorem 7.1 of [SS]). Therefore, 
it is enough to show boundedness for the sup norm of the zeroth, the first and the second derivatives of 
the functions on compact subinterval of (0, oo) for time coordinate t. 

Let us define the energy- momentum tensor T^p associated with the Lagrangian density (33), which 
has the form 

:= hABidaU^d^U'' - ^Vai^d^U^d^U''). (39) 

Components of T^/j are as follows, with (,) and || • || denoting the inner product and the norm with 
respect to h, 



= 





%0 


%^ 


Tet 












£ 


T 




T 


£ 


» 1 





- 


-t^g / 



\dtUf + \\deUf 
2{dtU,deU) 




2{dtU,deU) 
\\d,Ur + \\deUr 








t\\\dtUr-\\deUr) 



Note that equations for A can be simplify by these quantities as follows: 

dtX = t£, 
deX = tT, 
djX-djX = g. 



(40) 



(41) 
(42) 
(43) 



By using light cone estimate [MV] and Christodoulou-Tahvildar-Zadeh's identity [CTZ], we have the 
following lemma: 



Lemma 2 (Lemma 2 of [NM03]) There is a positive constant C such that 



£<C 



t e (0,oo), 



(44) 



where C depends only on the initial data at t = to. 



□ 



lemma 2, wc have 

-2(P+Pr)/ 



the desired bounds on | dtP |, 



deP 



-2{P- 



dtPi 



Proof of Theorem 1: From 
I dePi I, I e-2(«'^'+^/)at$j 

I g-2(p+p,)(5^g^ ^ ^Q^^^ _ 1^ I g-2(p+p,)(5^Q^ ^ ^Q^^^ _ 1^ for all t G (0,oo). Once we 

have bounds on the first derivatives of P, Pj and it follows that P, Pj and 4> arc bound for all t G (0, oo). 
Then, we have bounds on dt^i, 9t*, S^*, dtQi + ^dt^i - $/St* and deQi + '^de^i - ^/Sg*. 

Consequently, $/, dtQi and dgQi are bounded. Finally, we have boundedness on Qj. 

Next, we must show bounds on the second and higher derivatives of the functions. There is a well- 
known general fact that, in order to ensure the continuation of a solution of a system of semi- linear wave 
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equations, it is enough to bound the first derivative pointwise. Then, we have boundedness of the higher 
derivatives. 

By the constraint equations (41) and (42), boundedness and compatibihty with the periodicity in 6 
for the function A is also shown. The arguments for that are the same with one of the proof of Theorem 
1 of [NM03]. Thus, we have completed the proof of Theorem 1. □ 

4 Global existence theorem in constant mean curvature (CMC) 
time coordinate 

To show a global existence theorem in CMC time coordinate for T^~^-symmetric spacetimes, Hcnkel's ob- 
servation is applicable [HO] . He has considered prescribed mean curvature foliations in locally ?7 ( 1 ) x ?7 ( 1 )- 
symmetric four-dimensional spacetimes with matter fields. The system of our D-dimensional spacetimes 
is similar with one treated in that paper^. One difference between Henkel's and ours is target spaces of 
the wave maps. Fortunately, arguments in the proofs of Propositions 5.4-5.6 of [HO] do not depend on 
properties of the target spaces. In addition, we need not estimate functions of matter fields separately 
since the field equations of the matter fields are included in the wave map system in our case. Thus, from 
negativity of mean curvature of areal time slices. 



K=-e 



P-A 



< i^<0' ie(0,oo), (45) 



(where equation (41) has used) the spacetimes has crushing singularity into the past and thus a neigh- 
borhood near the singularity can be foliated by compact CMC hypersurfaces [GC]. Once we have one 

compact CMC hypcrsurface with negative mean curvature, it is shown that the CMC foliation covers the 
entire future of the initial CMC hypcrsurface (see [ARR]). Thus, the following theorem is obtained: 

Theorem 2 Let {M.D,gij,u, A.^, (p) he the maximal Cauchy development of smooth T^~'^ -symmetric Cauchy 
data on M.d-i ~ T^~^ for the T^~'^ -symmetric EMD system. Then, under the assumption (37), A4d 
can be covered by hypersurfaces of the constant mean curvature in the range (— oo,0). □ 



5 Existence of asymptotically velocity-terms dominated (AVTD) 
solutions 

5.1 Analytic case 
5.1.1 Generalities 

As a method to construct AVTD singular solutions, the Fuchsian algorithm has been developed [KR]. To 
answer the question rigorously whether the AVTD behavior exists or not for vacuum Einstein gravity in 
arbitrary spacetime dimension, the algorithm will be used. 

Now, we briefly review the Fuchsian algorithm. Let us consider a hyperbolic PDE system, 

F[u{t,x")]=0. (46) 

Generically, u can have any number of components. Here, we will assume that the PDE is singular 
with respect to the argument t. The Fuchsian algorithm consists of three steps: At first, identify the 
leading (singular) terms Mo(t, x) which are parts of the desired expansion for u. This means that the most 
singular terms cancel each other when uo{t,x) is substituted in equation (46). One can get the solution 
uo from a system of velocity-terms-dominated (VTD) equations which are obtained by neglecting spatial 
derivative terms from full equations (46). Second, introduce a renormalized unknown function v{t,x), 
which is given by 

u = Uo + t'^u. (47) 
^Rather, our system is simpler than Henkel's because there is no shift vector in ours. 
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If Mo ~ t*^, we should set to = A: + £, where e > 0. Thus, m is a regular part of the desired expansion for 
u. Finally, obtain a Fuchsian system for u by substituting equation (47) in equation (46). That is, 

{V + N{x))u = t"V{t,x,u,d^u), (48) 

where T) := tdt and TV is a matrix which is independent of t and a > 0. Note that one can always take 
a = 1 by introducing as a new time variable. V can be assumed to be analytic in all of arguments 
except t and continuous in t since the following existence theorem (Theorem 3) is a singular version of 
the Cauchy-Kowalewskaya theorem essentially. Note that equation (48) is a singular PDE system for the 
regular function v. 

Once we have the Fuchsian system, we can show the existence of a unique solution for prescribed 
singular part Uq by the following theorem. 

Theorem 3 (Theorem 3 of [KR]) Let us consider a system (48), where N is an analytic matrix near 
X = 0, such that \\(J^\\ < C for < (t < 1 (boundedness condition) and V is analytic in space x and 
continuous in time t. Then the Fuchsian system (48) has a unique solution which is defined near a; = 
and t = 0, and which is analytic in space x and continuous in time t, and tend to zero as t ^ 0. □ 

Note that the boundedness condition holds if every eigenvalue of A is non-negative. Theorem 3 implies 
that renormalizcd unknown functions must vanish as t ^ if the conditions are satisfied. Therefore, 
the only singular terms which are solutions to VTD equations remain and they are solutions to full field 
equations at t = 0. Thus, we obtain AVTD singular solutions. 

5.1.2 Application to vacuum T^^-symmetric spacetimes 

As mentioned in Section 1, it is interesting for us in the vacuum case. Let us take $/ = \l/ = (/) = a = 0. 
Then, we obtain a metric of the target space of the wave map from (38), 

D-l , D-l 



h = dP^+Y^ dPf + Yl e'^'^'dQl (49) 

1=3 1=3 

By Lemma 1 (or equations (17) and (16)^) we have a system of evolution equations for vacuum Einstein 
gravity as follows: 

V^Pj-t'dlPj = -le-^iP+P')[{pQjf-{tdeQi)% (50) 
V'^Qi-t'dlQi = 2[V{P + Pi)VQi-i'{deP + dePi)d0Qi]. (51) 
Neglecting spatial derivative terms from equations (50)-(51), one can obtain VTD equations as follows: 

V^Pj = -^e-^(P+P')iVQjf, (52) 

V^Qi = 2V{P + Pi)VQi. (53) 

From these equations (52)-(53) we can find VTD solutions as follows: 

Pr^{t,e) = pio{e)\nt+pn{9), (54) 
Q7™(i,^) = qio{e)+t'^P''+P'°\n{6), (55) 

where po '■= J2iPi^- Hereafter we put kj := po + pio- Thus, the following formal solutions which have 
the leading terms P/"^^ and QY^^ are obtained: 

Piit,e) = py^^ + t''TTi{t,e), (56) 

Qi{t,0) = Qr'' + t^'"^i{t,e), (57) 

where e/ > and fc/ > 0. 



*In this case, the dilaton and the Maxwell equations (12), (19), (20) are automatically satisfied. 
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Let us be tt/ := (ttq/, ttu, 7r2/) := (tt/, Ptt/, f9e7r/) and R := {kqi, ku, K21) ■= (Ki^VKijidgKi). From 
these, wc define v := (7f/,K/)^. We have a Fuchsian system for v by inserting solutions (56) and (57) 
into the system (50) and (51): 

(V + N)v = fV{t,e,v,d0v), (58) 
where S = min{2A;/ — e/, 2(1 — fc/) — e/, e/} for any /, ^ is a vector-valued regular function and N is 



N 



k3 



(59) 



where 



iV„ := 







-1 



ej 2ei 




and Nkj := 



0-10 
2ki 




(60) 



which are independent of t. It is easy to see that every eigenvalue of N is non-negative and regularity of 
the right-hand-side of the system (58) holds if and only if < e/ < min{2A;7,2(l — fc/)}, i.e. < fc/ < 1. 
Thus, we have the following theorem from Theorem 3. 

Theorem 4 Let pia{0), qiAifi), where A = Q, 1, be real analytic functions of 9 in a neighborhood of 
6 = 60 and < kj < 1 for all 0. Then, for the vacuum Einstein equations (50) and (51) in arbitrary 
dimensions (D > A), there exists a unique solution with the form (56) and (57) in a neighborhood of 
= 9o, where tt/, kj tend to zero as t 0. □ 

Theorem 4 implies that VTD solutions P^^^ {t, 0) and Q^'^^ {t, 0) are solutions to full Einstein equations 
(50) and (51) near the initial singularity. The VTD solutions depend on the maximal number (i.e. 
4 X (Z) — 3)) of arbitrary functions pia{9) and qjA{0)- Thus, it has been shown that T'^~^-symmetric 
vacuum spacetimes on T^~^ x R have AVTD behavior in general. 



5.2 Smooth case 
5.2.1 Generalities 

The argument in this subsection has been developed by Rendall [RA] (which was generalized in [SF]). 
Wc will review it, which is a smooth version of the Fuchsian algorithm. 

We want to extend to the C°° category of the existence of AVTD singular solutions. To do this, we 
introduce notions of regularity and formal solutions of equation (48). 

Definition 1 A function f{t,x) from an open subset fl C [0,00) x R" to R™ is called regular if it is 
C°° for all t > and if its partial derivatives of any order with respect to x gK"" extend continuously to 
t — (within Cl). 

Definition 2 A finite sequence (vi,. . . ,Vp) of functions defined on an open subset il C [0,oo) x R" is 
called a formal solution of order p of the differential equation (48) on Q. if 

1. each Vi is regular and 

2. {V + N[x))vi — tV{t,x,Vi,dxVi) = 0{f) for all i as t ^ in ^l, where O-symbol is taken in the 
sense of uniform convergence on compact subsets. 

If functions N and V in (48) are analytic ones, existence of analytic solutions can be shown [KR]. 
Contrary, if TV and V are merely smooth, the following existence theorem for formal solutions of any 
order is obtained: 
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Lemma 3 (Lemma 2.1 of [RA]) If V is regular and N is smooth and satisfies \\(J^\\ < C for some 
constant C and for < a < 1 in a neighborhood of zero, then equation (48) has a formal solution of any 
given order which vanishes at t = 0. □ 

The Fuchsian system (58) satisfies the conditions in Lemma 3 if e/ < 2fc/ and e/ < 2 — 2ki. Then, this 
system has a formal solution of any order which vanishes at t = 0. 

Now, we will rewrite the system (48) to a symmetric hyperbolic system defined as below: 

Definition 3 We say that the system of differential equations 

{A'^{t,x)V + J\f{x) + tA^{t,x,v)dj}v = tV{t,x,v) (61) 

is regular symmetric hyperbolic (RSH) if is uniformly positive definite and symmetric, the are 
symmetric, and all coefficients are assumed to be regular. 

For generic symmetric hyperbolic systems, it is a key point that energy inequalities defined in the L^- 
Sobolev space hold. The basic idea used to prove local existence results for symmetric hyperbolic systems 
is as follows. A family of approximated systems is constructed and is solved, and then, energy estimates 
are used to show that solutions of the approximated systems converge to a solution of the original systems 
in a certain limit [FR, MS, TM] . According to this, a sequence of analytic data as smooth and approximate 
data was used and convergence of the corresponding analytic solutions to a smooth solution was shown 
for the RSH system (61) coming from the Fuchsian system (48) in [RA]. 

It is impossible to get a tractable RSH system from the Fuchsian system in any cases, since when 
the original Fuchsian system would be rewritten to a symmetric hyperbolic system, two disadvantages 
may happen. One is violation of the boundedness condition for J\f (positivity of eigenvalues of Af) and 
another is change of powers of t in the right-hand-side of the system. To overcome the first problem, we 
will use the notion of formal solutions (and Lemma 3) and, for the second problem, we will consider as 
the cases may be. 

A formal solution of order p of the RSH system is defined as Definition 2. Note that it is necessary 
to verify that a formal solution of (58) of order p is also a formal solution of (61) of order p. If this can 
be done, we can proceed the next step. 

Given a formal solution {vi, . . . ,Vi} of (48), it is possible to consider the difference between a genuine 
solution V (which has been known yet) and the formal solution. Put Zi := t^~'^{v — Vi). Then, we can 
obtain a system for Zi as follows: 

{A°{t,x)V + [N'{x) + H-l)A°{t,x)\+tA^{t,x,Wi + t'-'^Zi)dj}zi = tVi{t,x,Zi), (62) 

for some regular function Vi. Note that choosing i large enough makes real part of the eigenvalues of 
J\f + {i — 1)A° positive, since A° is uniformly positive definite. 

If the data are approximated by a sequence of analytic data Sm '■= {piAm,QiAm), a corresponding 
sequence of analytic solution is obtained. We can construct the sequence Sm which converges to smooth 
data 5* := (piAjQia) in , uniformly on compact subsets. If the (approximate) formal solutions are 
constructed as in the proof of Lemma 3, then v^i ^ Vi as m ^ oo, uniformly on compact subsets. Here 
Vmi are analytic formal solutions of order i corresponding to the analytic data Sm and Vi is a formal 
solution order i for the smooth data S. The same is true for the spatial derivatives of these functions of 
any order. Therefore, one can conclude that the sequence of coefScients of the system (62) also converge 
on compact subsets as m — > oo. 

Now, our task will be show that the sequence of solutions Vmi exists for a time interval independent 
of m and has a limit as m — > oo solving the system (61). Fortunately, the global existence theorem 
(Theorem 1) implies that there is a sequence of smooth solutions to (62) on a common time interval for 
all m. Then, the following theorem ensures existence of genuine smooth solutions. 

Theorem 5 (Section 4 of [RA] and Theorem 5.1 of [SF] (see also Theorem 2.3 of [KS])) Let 

Zm{t,x) be a sequence of regular solutions on [0,ti) xU C [0, oo) x R", with Zm{0,x) = 0, to a sequence 
of RSH equations 

{A°m{t,x)V + Nmit,x)+tAi,{t,X,Zm)dj}Zm=tVm{t,X,Zm)- (63) 
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Suppose that is positive definite for each m and that the coefficients converge uniformly to A'^, Nqo, 
Ai^ and Voo on compact subsets as m —> oo, with the same properties as Am, N^, A^^ ^m, o,nd that 
the corresponding spatial derivatives converge uniformly as well. Then Zm converges to a regular solution 
Zoo of the corresponding system with coefficients Al^, Nqo, ^'^^ ^oo on [0,to) x for some to, and 

Here, we have fixed a value of i large enough and omitted the index i. 

Finally, we must verify that the obtained smooth solution of (61) is also a solution to (58). When the 
above arguments be done, the proof of the existence of smooth VTD solutions completes. 

5.2.2 Low velocity case 

Now, we will rewrite the Fuchsian system (58) to a RSH system as follows: 

{V + J\f + tA''dg)v = t<V{t, 9, V, dev), (64) 

where C = min{2fc/ — cj, 2(1 — ki) — ej, 1 — 2ki + ej, ej} for any I and J, V is a vector-valued regular 
functions and 



(65) 



where 








-1 




' 


-1 






2ei 


and Nk, ■= 





2ki 







-1 







-1 



(66) 



and 



A' = 



iHer 



A 



Her 



where ^Her 





0-1 
0-10 



(67) 



which are independent of t. In our case, A'^ becomes the identity. The regularity of the right-hand-side 
of (64) is satisfied if < fc/ < 3/4 for any /. 

We must verify that solutions to the Fuchsian system (58) are also solutions to the RSH system (64). 
The system (64) is obtained by replacing tdoTTai and tdgHQi with n2i and K27, respectively. For analytic 
solutions, the system (58) implies VitdgiTQi — 1:21) = and ^^{tdeKoi — K21) = 0. Since tdgnoi — 7^21 = 
and tdeKoi — K21 = at t = 0, we can conclude that these equations are satisfied in all time. Thus, an 
analytic solution of (58) is an analytic solution of (64). Repeating this argument for formal solutions of 
order i to the Fuchsian system (58), we have the same conclusion for formal solutions. 

Next, we want to apply Theorem 5 to our system (64). As we can see from (62) and (66), it is enough 
to take i > 3, where i is the order for the formal solution. Then, we have a smooth solution to the RSH 
system (64). 

Finally, we have to show that the obtained solution to (64) is also a solution to the Fuchsian system 
(58). From (64) and the fact that n and R arc vanish at t = 0, wc have 7121 — tdgTioj = c]{6)t and 
K21 — tdgKoi = Cj{0)t for some functions c]{0) and Cj{6). Suppose analytic functions converging to 7f 
and K be Wm and Then, we have 

sup I Cj{6)t I = sup I {1721 - tdgwoi) - (7I"2/m - tdgTTQim) \ 

e e 

< sup I TT2I - T^2Im | + SUp i | de^Kol - deTTQIm H 0, 
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for any t asm ^ oo. Here, equations TT2im — tdeTToim = for any m has been used. Thus, tt2I — tdenoi = 
was shown for all time. K2i — tdeKoi = can be also proved similarly. Thus, it was shown that the smooth 
solution of the RSH system (64) is also a smooth solution to the original Fuchsian system (58). To sum 
up, 

Theorem 6 Suppose that piA and qia o,re smooth functions of and < /c/ < 3/4 for all 9. Then, for 
each spatial point 9 = 9a, there exists a solution of the Einstein equations (50) and (51) in a neighborhood 
of 9 = Oq of the form (56) and (51) where 2ki — 1 < ej < min{2A;/, 2 — 2ki} for any I and J, and nj and 
Ki are regular and tend to as f — > 0. Given the form of the expansion and a choice of a, the solution 
is unique. □ 



5.2.3 Intermediate velocity case 



We want to complement the whole range < fcj < 1 arising in Theorem 4. Now, the following candidates 
for solutions to the vacuum Einstein equations will be considered: 



Pf (t, 9) = P7™ + a/(0)i2-2fe/ + t^-^'"+"7ri{t, 9), 



where e/ > and fc/ > and ai are chosen as follows: 

= I — I exp 



2 \2-2ki 



D-l 



2(p/i + ^Pji) 



7=3 



(69) 
(70) 



(71) 



Therefore, the number of arbitrary functions is 4 x (Z) — 3) still. Solutions (70) are the same with (57). 
Substituting (69) and (70) into the vacuum Einstein equations (50) and (51), the following RSH system 
is obtained: 



{V + N + tA^de)v = t^Vit, 9, v, dgv), 



(72) 



where ( — min{2 — 2ki — ej,3 — 2(fc/ + kj) + ej, 2ki — 1 — ej} for any / and J, V is a vector-valued 
regular functions and 



(73) 



where 



-10 
(2 - 2ki + tif 2(2 - 2ki + £/) 
0-1 



(74) 



which are independent of t. From > 0, we have a regularity condition 1/2 < kj < 5/6. Note that the 
choice for aj (71) is to eliminate the leading order terms in the RSH system (72) obtained after using 
the solutions (69) and (70). 

The procedure of replacement {tdenoi 7^21 and tdeKoi K21) is the same with one of the low 
velocity case (Subsection 5.2.2). Thus, we have the following result: 

Theorem 7 Suppose that pja and qjA are smooth functions of 9 and 1/2 < /c/ < 5/6 for all 9. Then, for 
each spatial point 9 — 9q, there exists a solution of the Einstein equations (50) and (51) in a neighborhood 
of 9 — 9o of the form ( 69) and (10) where 2(ki + kj) — 3 < ej < min{2A;/ — 1, 2 — 2ki} for any I and J, 
and TT/ and kj are regular and tend to as t — > 0. Given the form of the expansion and a choice of ei, 
the solution is unique. □ 
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Remark 1 It is impossible to cover the whole range < kj < 1 in the analytic case (Theorem 4) by 
Theorem 6 and Theorem 7, which are of smooth cases. Fortunately, we can overcome this problem by 
repeating the above method n times. The argument is the same with Section 5.7 of [SF], so the details 
will be omitted. If solutions (69) are replaced by 

n 

the regularity condition of the right-hand-side of the symmetric hyperbolic system becomes 

1 - 2[{n + 1) - {ki + nkj)] < ej < min{l - 2n(l - ki), 2 - 2fc/}, (76) 

for any I and J. Here, each ajj is defined as the leading order terms are canceled at each stage like (71). 
Therefore, the number of free functions remain 4 x (£) — 3) still. Then, we have the following inequality 
for ki: 

1- 7^ < fc/ < 1- ^ (77) 
2n 2(n + 2) ^ ^ 

Thus, the range (1/2,1) is covered by the infinite sequence of intervals (1 — l/2n, 1 — l/(2(n + 2))). 
By combining this result with Theorem 6, an existence theorem of smooth AVTD solutions under the 
condition < kj < 1 is shown. □ 



6 Conclusion 

We would like to comment on the structure of the spacetime into the future direction since we have not 
discussed on it. Prom point of view of the SCC, we want to show future completeness of any causal 
geodesic of the spacetime. Recently, it has been shown that Gowdy and ?7(l)-symmctric (in the case 
of small initial data) spacetimes are geodesically future complete [CB, CBC, CBMOl, RH]. One of key 
ingredients is to show energy decay by using corrected energy method. Fortunately, estimate of energy 
decay for our wave map can be shown [NM04] . Another ingredient is the geometric stnicture of target 
spaces. Unfortunately, our understanding of the structure of (35) is entirely out of reach at the present 
time. 
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